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ABSTRACT 



We analyze the rise and fall times of type la supernova (SN la) light curves discovered 
by the Sloan Digital Sky Survey-II (SDSS-II) Supernova Survey. From a set of 391 light 
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curves /c-corrected to the rest-frame B and V bands, we find a smaller dispersion in the 
rising portion of the light curve compared to the decline. This is in qualitative agreement 
with computer models which predict that variations in radioactive nickel yield have less 
impact on the rise than on the spread of the decline rates. The differences we find in the 
rise and fall properties suggest that a single 'stretch' correction to the light curve phase 
does not properly model the range of SN la light curve shapes. We select a subset of 105 
light curves well observed in both rise and fall portions of the light curves and develop a 
'2-stretch' fit algorithm which estimates the rise and fall times independently. We find 
the average time from explosion to S-band peak brightness is 17.38 ± 0.17 days, but 
with a spread of rise times which range from 13 days to 23 days. Our average rise time 
is shorter than the 19.5 days found in previous studies; this reflects both the different 
light curve template used and the application of the 2-stretch algorithm. The SDSS- 
II supernova set and the local SNe la with well-observed early light curves show no 
significant differences in their average rise-time properties. We find that slow-declining 
events tend to have fast rise times, but that the distribution of rise minus fall time 
is broad and single peaked. This dist ribution i s in co ntrast to the bimodality in this 



parameter that was first suggested by IStrovinkI (120071 ) from an analysis of a small set 



of local SNe la. We divide the SDSS-II sample in half based on the rise minus fall 
value, tr — tf < 2 days and tr — tj > 2 days, to search for differences in their host 
galaxy properties and Hubble residuals; we find no difference in host galaxy properties 
or Hubble residuals in our sample. 

Subject headings: supernovae: general, methods: data analysis 



1. Introduction 



Type la supernovae (SNe la) are bright stellar explosions important for their role a s distance 



indicators. SN la distances have been used to constrain the value of the Hubble constant ( Jha et al 



1998: Freedman et al 



200lh. and they also showed that o u r universe has a lo wer than critical mat 



ter de nsity (jGarnavich et al.lll998l : |Perlmutter et al.lll998l ). iRiess et all (|l998l ) and lPerlmutter et al 



(| 19991 ) used distant la SNe to show that the universe currently has an accelerating rate of expan- 
sion implying a 'dark' energy component. As more SN la observations h ave confirmed this result 



(|Knop et al.ll2003l : lTonrv et al 



2003 



Harris et al 



constraining the properties of dark energv (iGarnavich et al 



20071 : IWood-Vasev et al 



2007 



20061: 



Miknaitis et al 



2007 



Riess et al.l 12004'). the focus has shifted to 
1998bl: kstier et al. 



Eisenstein et al 



2OO7I ). Systematic uncertain- 



20061 : 



Riess et al 



ties in measuring da rk energy parameters with supernovae now dominate over statistical errors 
(jKessler et al.ll2009al ) and a better understanding of supernova physics may help to constrain dark 
energy properties using SNe la. 



It is widely accepted that SNe la are the result of thermonuclear explosions of carbon-oxygen 
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white dwarfs (WDs), but the nature of the progenitor remains uncertain. In most models, the 
explosion occurs when the WD nears the Chandrasekhar limit by gaining mass from a binary 
companion. This mass gain is achieved either by single-degenerate (SD) mass transfer or WD 
coa lescen c e in a double-degenerate (DD) scenario. For a comprehensive review of these models, 



see iLivid (l200d) . Observers have attempted to distinguish between the two models, with con- 
flicting results. iHowell et al.l (|2006l ) found that the extremely luminous supernova SNLS-03D3bb 
(SN 2003fg) had low ejecta velocity that could have resulted from a super-Chandrasekhar progen- 
itor. The high total mass and larg e nickel yield suggested that SNLS-03D3bb was the product of 
a DD merger. iHicken et al.l (|2007l ) found that SN 2006gz had attributes consistent with the DD 
model. Its spectrum showed significant amounts of unburned carbon, and a low silicon velocity at 
early phases. T he very broad ligh t curve implies a large yield of radioactive nickel as expected in 
the D P model (IHicken et al.l 120071 ) . although late-time observations show only weak iron emission 
lines (iMaeda et al.ll2009l ). Even though the SD model is generally accepted as the most plausi- 



2007 



ble (iLivid l2000l) , there is evidenc e that it might not be the complete story ([Hicken et al. 
Pritchet. Howell. &: Sullivanll2008l ). Some mixture of progenitors may be producing events with 



subtle differences in character. The explosion rates for these progenitor systems will depend on the 



tors may not be constant through the history of the universe 


(Hamuv et al. 


1996bl: 


Gallaeher et al. 


2005; 


Sullivan et al. 


2006: 


Mannucci. Delia Valle. h Panaeia 


200d: 


Gallagher et al. 


20081). 



Understanding t he origin of th e Phillips relation (the co rrelation between light curve shape 
and peak luminosity; jPhillipd (|l993l ): iKasen &: WooslevI (120071)) and its dispersion is critical to im- 
proving SNe la as reliable distance indicators. iTimmes et al.l (|2004l ) argued that radioactive nickel 
yield should be determined by progenit or metallicity, but age of the host galaxy stellar popula- 
tion appears to control the mass of ^^Ni (IHamuv et al.lll996bl: lOallagher et al.ll200,4 ISullivan et al 



2006). Progenitor metallicity m ay play a secondar y role ( Gallagher et al. 20081') or ha,ve no s ignifi- 



cant effect on radioactive yield (jHowell et al.l l2009l ) . Recent models by IWooslev et al.l (|2007l ) show 
that variations in kinetic energy (KE) , metallicity, and mixing between burning layers provide light 
curves beyond the range of the observed Phillips relation, implying that not all combinations of 
variables are found in real SNe la. 

SN la luminosity and light curve shape may be influenced by the physics of the explosions. 
There is a consensus that normal SNe la result from a detonation (supersonic burning) of much 
of the progenitor WD, but it likely begins as a deflagration (subsonic fusion f ront) to al l ow ex - 
pansion and some bu r ning at low densities. This model was first introduced bv iKhokhlovl (jl99ll ); 
Hoefiich &: Khokhlovl (|l996l ) showed that the model matched the observed light curves and spectra 
very well. In the "delayed detonation" scenario, the timing of the transition from deflagration 
to detonation provides a natural means of varyin g the radioactive nick el yield and generating the 
observed range of luminosities and decline rates (jArnett &: Livndll994l ). Deflagrations tend to be 
very asymmetric, and if the asymmetries surviv e the detonation then viewing angle can result in 
perceived variations even for similar explosions (jKasen &: Plewall2007l ). 
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Study of the early -time light curve may be im portant in dia gnosing the proge nitor problem 
and explosion physics (jHoefiich fc Khokhlovl Il996l ) . The work by iRiess et al.l (|l999l ) to constrain 
the early light curves of local supernovae established the decline-rate-corrected average rise time of 
19.5 ±0.2 days. Her e , "ris e tim e" is defined as the time elapsed from explosion to peak i?-band 



flux. lAldering et al.l (|200d ) and iGoldhaber et al.l (|200ll ) demonstrat ed the consis t ency between 



high- and low-reds hift rise times, and similar results were found bv IConlev et al.l ()2006l ) with a 
large set of SNLS (jPritchet et al.ll2005l ) supernovae. lOarg et al.l jSoJ' studied the rise times of 
SNe la discovered behind the Large Magellanic Cloud during the SuperMACHO survey lens project 
and measured a 17. 6±1. 3-day rise time in the equivalent of the V bandpass. 



Recently, IStrovinkI (|2007l ) analyzed eight low-redshift SNe la with well-observed early light 
curves from a new perspective. A single "stretch" parameter has commonly been used to describe 
the full range of B and V light curves by coni pr essing or expan ding the time axis around the 
epoch of peak brightness (jPerlmutter et all 119971 ). IStrovinkI (120071 ) decoupled the rise and decline 
portions of the eight nearby light curves and found that SNe la can have a range of rise times for 
a given decline rate. His rise-time minus fall-time distribution (t^ — tf, hereafter defined as RMF) 
roughly divided his small supernova sample into two groups, possibly suggesting two progenitors 
or explosion mechanisms. A larger set of SNe is required to rigorously test this possibility. 

Here, we arialyze light curves from the Sloan Digital Sky Survey-II (SDSS-II) Supernova Survey 
(jFrieman et al.ll2008l ) which spectroscopically identified approximately 500 SNe la over its three- 
year lifetime. The SDSS-H Supernova Survey scanned 300 deg^ of sky with a cadence as rapid as 
2 days between visits (weather and lunar phase often increased the time between observations of 
the same field) , making the survey well suited to an early rise-time study. We also introduce a new 
fitting method that independently estimates the rise and the fall times of SNe la light curves. 



Data 



2.1. SDSS-II Light Curves 



The SDSS-II Supernova Survey (jFrieman et al.ll2008l ) was designed to find and characterize 
several hundred SNe la at intermediate redshifts in order to fill in the supernova "desert" between 
the nearby discoveries and the"high-2;" events. Mapping the expansion his tory at z ~ 0.2 pro- 



vides unique tests of cosmological rnodels and constrains systematic errors (jKessler et al 



2009a 



Sollerman et al.ll2009l : lLampeitl et al.ll2010l ). The SDSS-II supernovae also provide a large, uniform, 
high-quality sample of SNe la to study the properties of these explosions. The SDSS-II Supernova 
Survey operated thr ee campaigns betw een 2005, 2006 and 2007 September and December using 
the 2.5 m telescope (iGunn et al.l 120061 ) at Apache Point Observatory to scan 300 deg2 of sky as 
often as every-second night. Template images were subtracted from each new night of data and 
software scanned for new variable objects. Every candidate supernova was inspected visually to 



avoid image artifacts and asteroids. Transients with a high probability of being SNe la (jSako et al 
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20081 ) were queued for spectrographic observation (|Zheng et alj 120081 ) . The SDSS-II SN survey 
takes advantage of the extensive data base of refe r ence images, object catalogs, and photometric 
cahbration compiled by the SDSS (see lYork et alJ (120001) for an overview of the SDSS). The early 
data release of the SDSS can be found inlStoughton et alJ (20021 ). Data release seven, the final data 



release of SDSS-II, can be found in lAbazajian et alJ (|2009l ) 



The superi iova sample miniin izes photometric errors by using a single f ilter set and a photo- 



1996 



metric system (ISmith et alJ |2002| ) calibrated by the origi nal SDSS project (jFukugita et al.. 

Gunn et al. 19981 ). The phot ometric qua l ity is assessed in Ivezic et al. ( 2004 ). and the astrometric 
calibrations are described in iPier et al.l (l2003l). Specifically, t he area of the SDSS-II SN survey 
has been calibrated to 1% photometric error (jSmith et al.ll2002l ). The photometric calibration em- 
ployed severa l mino r telescopes as well; this Monitor Telescope Pipeline (MTPIPE) is described in 
Tucker et al] kood ). 



A software "robot" was developed that automatically reduced CCD observations of standard 
stars in orde r to assess photo metricity and build data on site conditions at the Apache Point 
Observatory (jHogg et al.ll200ll ). Using this tool to select the best observing night s, weather and 



lunar phase extend the time between observations to an average of about 4.5 days (jFrieman et al 



20081 ). The rapid cadence of the SDSS-II SN survey assures that most of the SDSS SNe have 
well-defined, densely sampled light curves even on the rise portion. The data include pre-explosion 
flux measurements, so that many of the SDSS-II Supernova Survey type la light curves are among 
the most well-sampled and well-calibrated light curves that have been studied to date. 

The SDSS-II SN survey identified 498 spectroscopically confirmed SNe la (number accurate to 
the time of our analysis). Five o f these are extre mely peculiar and are not included in this study. SN 
2007qd (|McClelland et aklboiol ) and SN 2005hk dPhillips et aklboOTl ) are categorized as 2002cx-hke 
objects. SN 2002cx showed evidence for a very low expansion velocity compared to 'normal' SNe la 
as well as an unusually low peak luminosity , leading to the hypothesis that it was a pure defjagration 



event ([Branch et al 



2004 



Jha et al.ll200fil ). SN 2005gj, SN 7017, and SN 15557 (jPrieto et al.ll2007l ;i 



are 2002ic like. SN 2002ic showed a type la spectrum along with strong hydrogen emission lines 
(jHamuy et al.ll2003l ) suggesting a shock interacting with circumstellar material. 



We limit the events to redshift less than 0.4 and require that there be at least one photometric 
epoch more than 2 days before maximum and more than 5 days after maximum. These require- 
ments result i n a sample of 391 supe rnovae. The photometry is measured using "scene modeling" 
developed by Holtzman et al. (j2008l ). The Sloan q, r, and i photometry has been ^-corrected to 
Bessell B and V bands using SNANA (IKessler et al.ll2009bl : Hsiao et ahlboOTI ). We use the time of 
maximum obtained by the fitting procedure in SNANA as the initial value for our minimizing 
function. The redshift of the 391 SNe la ranges from 0.037 to 0.40 with a median redshift of 0.21. 
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Fig. 1. — : Top left: individual SDSS-II supernova observations plotted in the rest frame and normalized to a peak 

flux of unity. The light curves were fc-corrected to B band using SNANA. Ovcrplotted are the MLCS2k2 fiducial 
light curves with a 16.8-day rise time. Top right: median of all the data points binned by day are plotted (-I-) and the 
gray regions show the rms deviations about the median values. The white lines show the MLCS2k2 fiducial curves. 
Bottom: fluctuations about the median arc plotted along with template light curve residuals. This clearly shows 
the scatter in rise time is significantly less than in the decline portion of the light curve. The solid line shows the 
MLCS2k2 template with a 16.8-day rise time is an excellent match to the median of the 391 SNe la. The dashed 
line indicates that increasing the rise time by 0.7 days is clearly a poor fit to the median rise data while a decline of 
16.4 days is a good match. 



-7- 



Figure [J displays the entire set of 391 B and V light curves after estimating the time of 
maximum, normalizing the peak to unity, and correcting for time dilation. The data were binned 
in 1-day intervals and the median value calculated for each day, also shown in Figure [TJ The 
MLCS2k2 fiducials are an excellent match to the median light curves, with an extrapolation of 
16.8 days. The fading portion of the SDSS-II i?-band data is brighter than the fiducial, indicating 
that the supernovae discovered by SDSS-II are, on average, slower than typical low-redshift events. 
For the V band, the fiducial and the SDSS-II light curves match well. It is clear that the scatter 
about the median light curve is smaller during the rise than during the decline. The maximum 
root-mean-square (rms) dispersion for B band in 1-day bins is 0.09 (scaled flux units) on the rise 
and 0.13 on the decline. The average rms about the median on the rise is 0.069 while on the fall it 
is 0.091. This implies a smaller range of rise times than fall times and therefore that the full range 
of SN la light curve shapes cannot be well fitted by a single-stretch parameter. 



KasenI (|2010l ) has provided a theoretical prediction of the impact on the light curve of SN la 
due to the interaction of the WD explosion with a potential companion star. It is estimated that 
approximately 10% of SN la explosions that occur in the SD channel would have light curves that 
are affected by the shock interaction of the explosion with a binary companion. This effect is not 
immediately obvious in the SDSS-II from any analysis in this paper; however, a n in-depth sta.tistica l 



study is in progress to rigorously test the SDSS-II SN la sample for this effect (jHavden et al.ll2010l ). 



2.2. The Low-Redshift SN la Set 



In order to compare our fitter with the results in Strovink (j2007l ). we applied our 2-stretch 
fitting method to the same set of eight low- z SNe la studied by Strovink. These supernova e 



are SN 1990N (Lira et al 



SN 2001el (Krisciunas et al 



20031). SN 2 002bo (iBenetti et al 



998h. SN 1994D (IPatat et al.lll996l). S N 1998aq (iRiess et al.l hOQH) 



2004), SN 2003du 



Stanisheveta 



20071 ). SN 2004eo (IPastorello et al.ll2007bl ) . and SN 2005cf (jPastorello et al.ll2007al ) . IStrovinkI (l2007l ) 
combines light curve observations from different studies in his analysis of low-z supernovae, but this 
has the danger of confiicting calibrations and differences in photometric error estimates. Instead, 
we limit observations to a single published data set selected based on cadence and number of data 
points before maximum. Table [J compares the Strovink and the 2-stretch fitter results. We find 
very simi l ar va lues for these supernovae. The average difference between the rise time measured by 
StrovinkI (j2007l ) and our 2-stretch fitter is —0.23 ±0.58 days (standard deviation) and the difference 
between fall times is 0.15 it 0.47 days (standard deviation). 



3. 2-Stretch Fitting Method 



With a handful of well-observed events, IStrovinkI (j2007l ) found that the rise time of SNe la can 
vary for a fixed decay time. He developed a fitting method called Adaptive QUArtic Algorithm 
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(AQUAA) which employs quar tic sphnes with fuzzy knots to create smooth UBVRI Hght curves 
for individual supernovae. The IStrovinkI (|2007l ) fitter is computer intensive and does not work as 
well with noisy data; both are drawbacks when dealing with the large quantity and diverse quality 
of SDSS-II supernova data (see Section 2.1). 

For the SDSS-II data, we desired a fitting method which maintains the simplicity of fitting 
light curves with the single-stretch technique, while decoupling the rise and fall portions during 
the fitting process. Our '2-stretch' method attempts to accomplish this goal by stretching the 
pre-maximum portion separately from the post-maximum portion of the light curve. Our template 
curve is required to stay continuous at the joining point: peak brightness. Maximum light has 
a zero first derivative, making this a natural cutting point for using a 2-stretch method. A 
minimization is performed to determine the four parameters of the fit: rise stretch (s^), fall stretch 
(s/), time of maximum, and peak flux. It is not possible to completely remove the covariance 
between rise and fall times, because photometric errors and sampling gaps produce uncertainty in 
the time of maximum which mixes rise and fall times. Fitting simulated light curves with known 
rise and fall times allows us to estimate the correlation between rise and fall errors. 

Our implementation of th e fitting method uses B and V fiducial curves generated by MLCS2K2 
(|.Tha. Riess. fc KirshneiiboOTl l for a typical, A = supernova. The MLCS method has the advantage 
of creating a template from a large number of nearby supernovae. However, few of the supernovae 
used to train MLCS2k2 had light curve information earlier than 10 days before maximum light. 
For the region earlier than —10 days, we assume a simple expanding fireball mo del which implies 



that the fiux from the moment of explosion increases as the square of the time (iGoldhaber et al 



2001 



Riess et al.lll999l ). We then create a set of possible early-time light curves that all join with 
the MLCS2k2 fiducial at —10 days (see Figure [2]). This set ranges in rise time from 14 to 20 days 
in half-day intervals. The first derivative of the joined curves is not required to be continuous. In 
Section 4.3, we describe how we use the SDSS-II supernovae to determine which of these template 
curves is the best match to the observations. 

Historically, there are a number of ways to characterize the light curve shapes of SNe la. For 



Table 1. Comparison of 2-stretch Method versus IStrovinkI (j2007l ) 



SN Amis (inag) Amis (mag) trise (days) trise (days) tr ~ tf (days) tr —if (days) 

Strovink07 2-Stretch StrovinkO? 2-Stretch Strovink07 2-Stretch 



SN 1990N 
SN 1994D 
SN 1998aq 
SN 2001el 
SN 2002bo 
SN 2003du 
SN 2004eo 
SN 2005cf 



0.990 ± 0.034 
1.344 ± 0.021 
1.042 ± 0.021 
1.168 ± 0.021 
1.162 ± 0.031 
0.982 ± 0.021 
1.403 ± 0.037 
1.068 ± 0.021 



1.029 ± 0.017 
1.422 ± 0.023 
1.089 ± 0.013 
1.080 ± 0.010 
1.168 ± 0.031 
1.013 ± 0.010 
1.326 ± 0.010 
1.045 ± 0.008 



20.01 ± 0.46 
15.39 ± 0.47 
17.52 ± 0.58 
18.00 ± 0.56 
16.05 ± 0.37 
17.71 ± 0.35 
16.64 ± 0.44 
16.62 ± 0.25 



19.81 ± 0.25 
15.43 ± 0.23 
17.13 ± 0.22 
18.22 ± 0.18 
17.50 ± 0.19 
17.77 ± 0.17 
16.62 ± 0.13 
17.31 ± 0.10 



4.04 ± 0.74 
2.04 ± 0.53 
2.03 ± 0.64 
3.57 ± 0.63 
1.57 ± 0.53 
1.64 ± 0.44 
3.86 ± 0.53 
1.35± 0.36 



4.31 ± 0.26 
2.73 ± 0.23 

2.20 ± 0.22 

3.21 ± 0.18 
3.21 ± 0.21 
2.11 ± 0.17 
3.39 ± 0.13 
1.97 ± 0.10 
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Fig. 2. — : Comp arison of the MLCS2K2 fiducial curves (solid), with the Leibundgut template (dashed) used by 
Riess et all (|l999h . The difference between the Leibundgut template and our fiducial curve is 1.8 days at about 0.66 
mag on the risin g portion, w hich is important in explaining the approximately 2-day difference between our rise 



time and that of iRiess et al.l (|1999). Also, application of the single-stretch method using a very wide template will 
'squeeze' the light curve and push the maximum to a later date in order to fit the real rise data. The extrapolations 
shown here are described in Section 4.3; this figure displays 7 of our 13 extrapolations to the MLCS2k2 template. 



example, the MLCS2k2 -B-ba nd fiducial curve fall s 1.1 mag in 15 days after maximum, correspond- 
ing to a Ami5(i?) =1.1 mag (jHamuy et al.lll996bl ). Equivalently, the fall time, tj, is defined as the 
number of days it takes a light curve to fade 1.1 mag after maximum brightness. The 2-stretch 
method estimates fall stretch, Sf, and the two parameters are simply related hytf = Axsf, where 
^ is a constant that describes how long it takes the fiducial to fall 1.1 mag. The relation between 
Ami5(i?) and tf depends on the precise shape of the fiducial curve, but is approximately quadratic 
for MLCS2k2: 



tf = 14.711 - 9.631(Ami5(B) - 1.1) + 9.391(Ami5(5) - 1.1)^ 

over 0.7 <Ami5(i?)< 1.5 for the MLCS2k2 fiducial B band used here. We derived this equation 
by phase-stretching the MLCS2k2 fiducial light curve and directly measuring the tf and Ami5(i3). 
Note that we use the literal definition of Ami5(i?) and this may not be the same value derived 
by Ami5(i?) fitters (jPrieto et al.l l2006l ) precisely because rise and fall rates are not absolutely 
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correlated. 

In order to estimate the rise and fall times of each supernova, an routine was developed 
that modifies the fiducial curve, F{t), to fit the data at rest-frame times, t, by performing a 
minimization. Mathematically, the 2-stretch function S is defined as 



S{t) 



foF{{t-to)/sr) ift<to 
_ foF{{t-to)/sf ) if t>to 

, where Sr is the rise stretch, s/ is the fall stretch, to is the time of maximum light, and /o is the 
peak fiux. The fiducial curve is prepared so peak flux occurs at r = and the peak has a value of 
unity. The function S is an approximation to the true observed light curve and is then interpolated 
at each observation time. The minimization is performed on data out to 25 days past maximum 
light. 



We applied the 2-st retch fitter to the SNe la studied in iRiess et al.l (j 19991 ) and a subsample 
from iHicken et al.l (|2009l ). In all, we fit 41 i iearby SNe with the 2-stretch algo rithm; this sample 
includes the 8 SNe la analyzed bv IStrovinkI (j2007l ). 4 more SNe la a nalyzed bvlRiess et al. I (|l999) 
that are not in the IstrovinkI (12003) set, and 29 SNe la from CfA (jllicken et al.l bood ) that had 
pre-maximum data. A number of SNe are common to all three data sets, and in the se cases we 
use th e data with the highest number of pre-maximum observations. SN 1996bo from lRiess et al.l 
(| 19991 ) is poorly constrained on the rise portion of the light curve, and was not included in our 
analysis. Our minimizing func tion was unable to find a suitable minimum for SN 1997bq and 
SN 1998ef from lRiess et al.l (119991 ). Lastly, SN 2005hk was not used from the CfA set because it is 
a peculiar SN la, and has been omitted from the SDSS-II analysis as well. 



4. Method 



This section will explain issues not fundamental to the 2-stretch fitter itself, but are important 
issues post-fit in analyzing the data. First, we discuss our error calculation method, and then 
we explain how we use the errors to select our subsample of the 'best' SN la light curves. Next, 
we describe our method for template selection, ultimately resulting in the selection of a 16.5-day 
rise-time MLCS2k2 template which is used for all 2 stretch fits in this paper. Lastly, we discuss 
how we combine the B- and V- band fits for each SN la. 



4.1. Method of Estimating Errors 

The uncertainty in rise and fall measurements arises from both photometric errors and the 
distribution of observations across the light curve. To estimate the rise/fall error, we use a Monte 



Interactive Data Language, |http: / / www.ittvis.com/idl/idl7.asp | 
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Carlo method similar to that of IContardo. Leibundgut. &: Vaccal (|200d ). At each observation in a 



light curve, we employ a Gaussian distribution centered on the observed flux and with a standard 
deviation equal to the flux error of that point. We generate simulated photometry at each observed 
time and then perform the 2-stretch ^ minimization on the simulated light curve. The random- 
ization is performed 100 times, and the standard deviation of the set is used as the error for the 
rise/fall times. There is no distinguishable difference in error values if this is performed 100 times 
or 1000 times. The rise/fall- time uncertainties reported by this method are used as a measure of 
the quality of the fit of each SN. 

For the low-redshift supernovae, we found that the photometric uncertainties reported in the 
published light curves had a large range even though the supernova apparent peak brightnesses were 
similar. We assume this is due to some authors including systematic errors (photometric calibration 
uncertainties) in some analyses while other studies list only statistical photometric uncertainties. 
We add 0.01 in quadrature to all error values to take into account uncertainty in the model. 

We notice a bias toward larger errors on fast-declining supernovae. The reason for this is 
twofold: these supernovae are generally fainter for a given redshift and therefore have larger pho- 
tometric errors, and the rest-frame cadences of observations of these supernovae are longer. By 
excluding supernovae with large rise/fall error estimates, we have reduced the number of very fast 
decliners in our sample relative to supernovae with normal and slow-declining light curves. By 
a similar argument, we expect that the number of broad light curves is over represented in this 
sample when compared to the number of fast-declining events. 



4.2. Selecting Well-Sampled Light Curves 

To determine how the rise and fall times of individual SNe la are related, we had to choose 
supernovae with well-defined light curves before and after maximum light. Sorting supernovae 
based on the number of points within a certain time range does not take into account the variation 
in photometric quality at various redshifts. We decided to select light curves based on the rise- and 
fall-time errors calculated by our Monte Carlo method. 



Since one of our goals was to look for the double-peaked RMF distribution found by IStrovink 



(|2007l ) which has a separation of about 3 days, we decided to use a 2.0-day error cut on the rise time 
and the fall time. We applied these cuts to both the B and V bands before averaging the stretches 
from each band. These criteria produce a sample of 105 supernovae, out of the 391 SDSS-II SNe. 
The redshifts for these events range between 0.037 and 0.230; the average redshift is 0.14. 

Just as in the single-stretch method, it is possible to "correct" the observed light curve to 
best match the fiducial curve by dividing the time axis by the estimated stretch parameter. This 
normalizes the light curve to match the template curve that was used in the fitting process. The 
only difference for the 2-stretch method is that the pre-maximum and post-maximum time axes 
have different stretches applied. Figure [3] shows the 2-stretch-corrected light curves for the 105 
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SDSS-II light curves well sampled in both rise and fall. 




-30 -20 -10 10 20 30 
Phose (doy) 



Fig. 3. — : 105 SNe used in this study, corrected by both the rise and fall stretches, along with the flux and time 
of maximum found by the fitter using a 16.5-day rise-time template. The rise time of the best template is not the 
average rise time of the data, but a measure of the best relative shape between the extrapolated early rise portion 
and the MLCS2k2 template. This represents B band only. This normalized light curve has a jdoj of 211/247, a 
probability of fit of 0.988. This fit is much better than the single-stretch version, which has a probability of fit of 
0.242. 



4.3. Selecting the Fiducial Curve 

The MLCS2k2 fiducial light curves were poorly defined for early phases, so we constructed a 
set of extrapolated curves with a range of explosion dates that meet the MLCS2k2 curve at —10 
days. Recall that in the simple fireball model, the extrapolations are quadratic functions with zero 
flux at the day of explosion. We created extrapolations at half-day intervals starting with explosion 
at 20 days before B maximum and ending with explosions at 14 days before peak (see Figure [2|). 

To determine the best extrapolation, we fit each of the 105 light curves with the 2-stretch 
algorithm using each of the 13 extrapolated rise curves as a fiducial. We then divide the pre- and 
post-maximum portions of each light curve by the derived rise and fall stretches, and combine all 
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105 light curves into a single normalized light curve for each candidate fiducial curve. We calculate 
a value on the data points contained within the region of phase from —20 to —10 days. We 
found that the smallest reduced value and largest corresponding probability of fit were attained 
with the fiducial curve with a 16.5-day rise time in the B band (see Figured]). This template was 
best in V band as well. This MLCS2k2 16.5-day template is the template that was used for all 
analyses in this paper. We emphasize that this number is not our estimate for the typical rise time 
of a type la SN, but a determination of the shape of the early light curve that best matches the 
data, relative to the shape of the actual MLCS2k2 template data later than —10 days. We also 
applied a single-stretch fit and found that the probability of fit was always significantly smaller 
than for the 2-stretch fit, showing that the 2-stretch fitting method is a better model for SN la 
light curves (see Figure The value goes from 1515 (single stretch) to 1140 (2-stretch) with 
the addition of another model parameter, with 1050 and 951 degrees of freedom (dof), respectively. 
At the minimum x^j the equivalent probability of fit for the 2-stretch fitter is 0.988 while for the 
single-stretch fitter it is 0.242. 
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Fig. 4. — : Plot of the distribution for the fiducial curve rise times, consisting of 105 high-quality SDSS-II light 
curves. The light curves were normalized to the fiducial curve using the values determined by our 2-stretch fitter, 
then the reduced value of the region between —20 and —10 days was calculated. Plotting the reduced for 
the entire light curve fitting region produces a similar result in the comparison of single stretch vs. 2-stretch, but is 
slightly noisier. These data represent the B band only. 
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Lastly, we used our error calculations to determine the optimum day after maximum to cut 
off the fitting algorithm. The fraction of the 391 supernovae that pass our error cuts is a good 
estimator of the overall quality of the fits. We fit all of the light curves with a range of cutoff ages 
from +15 days to +35 days and found the best cutoff to be around +25 days. Longer age cutoffs 
can bias the fall time with noisier data that may not be as well represented by the template; we 
found that the best combination of small errors and negligible bias occurred at the cutoff age of 
+25 days. This age cutoff was used in all analyses reported in this paper. 



4.4. Combining the B- and F-Band Stretches 



Following IStrovinkI (|2007l ). the stretches from the B and V bands are combined to produce a 
single rise and a single fall-time estimate for each supernova. The stretches in the two bands {sr{B) 
and Sr{V), for example), found using the 16.5-day MLCS2k2 template to find the stretches, are 
combined using their weighted average, with our derived errors as the weights (in the form l/o"^). 
This has the advantage of reducing noise in the rise/fall-time measurements. Figure [1] shows that 
the SDSS-II B-hand light curves, on average, decline more slowly than the MLCS2k2 fiducial. 
This could be due to selection bias that comes from the tendency of a magnitude-limited search 
to find brighter than average events that also tend to have slower-declining light curves. However, 
the average SDSS-II F-band light curve is a good match to the MLCS2k2 IZ-band fiducial. This 
mismatch between the average B and V curves and the fiducial is puzzling. We have chosen 
to ignore this difference and simply renormalize the fiducial stretches to best match the SDSS-II 
average light curves in both bands. 

We found that the MLCS2K2 templates in B and V were different between their average rise 
and fall stretches when compared to the SDSS-II data. The y-band stretches on the fall portion 
were on average smaller than the B band, while on the rise portion the V-band stretches were 
larger than B band. In order to renormalize our template, we used the average stretches from the 
set of 105 SDSS-II SNe. We divide all B- and F-band stretches of the actual data by the average 
B-rise, i?-fall, V-vise, and V-fall stretches. This ensures that the templates have the same rise and 
fall stretches in both B and V, so that we can average the stretches together with confidence that 
they are centered about the same mean value. This method of modifying the stretches of the data 
after fitting relies on the fact that starting our minimizer at a different initial condition will 
still provide the same final result. The other choice is to modify the template before fitting, in 
which case we would be multiplying the template by these average stretches from our data set. As 
long as the minimizer finds the same minimum, these two methods are interchangeable. Testing 
indicated this to be the case. 

To clarify the relation between B and V light curves, the difference between the rise stretches 
and fall stretches after renormalizing the templates for the two bands is plotted in Figure [5l We find 
that on the rise, the B- and y-band stretch difference is constant with decline rate, suggesting that 
the shape of the rise portion of the fiducial curve is consistent between the two wavelengths. But the 
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plot of Sf{B) — Sf(y) shows a kink near Ami5(i?)=0.9, indicating that the fading portions of the 
B and V hght curves require two different stretch parameters to match the data. Specifically, very 
slow declining events tend to fade more slowly in the blue than in the visual band when compared 
to normal decliners. This effect should make very slow decliners more blue after maximum than 
their faster-fading cousins. 




0.5 1.0 1.5 2.0 

Am,5(B) 




0.5 1.0 1.5 2.0 

Am,5(B) 



Fig. 5. — : Left: the _B-band rise stretch minus IZ-band rise stretch for 105 SDSS-II type la SNe. The selection 
criterion for these 105 is described in Section 4.2. The B- and l/-band rise stretches show no clear variation across 
the sample, indicating that the rise stretch is consistent between the two bands. Right: the difference between the 
B-band fall stretch and the '[/-band fall stretch for the same 105 SNe. Here, the data show evidence for variation 
with Ami5(_B). Specifically, we find that slow-declining events have a larger than expected _B-stretch than V^-stretch 
when compared to events with Ami5(_B)> 1.1. 



Clearly, the i?-band fall stretch is not strictly correlated with the F-band fall stretch and this 
may indicate a physical difference between the slowest fading events and more normal-declining 
supernovae. However, the effect is small enough that we do average the derived stretches for this 
analysis. 



5. Analysis 

5.1. Rise versus Fall Times in the SDSS-II 

Using our sample of 105 high-quality SDSS-II supernovae, with the 16.5-day MLCS2k2 tem- 
plate for the fits, we obtain an RMF distribution that is best described by a single Gaussian 
distribution. This distribution can be seen in Figure O which shows a standard histogram of the 
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RMF distribution along with an ideogram of the distribution q These histograms, along with the 
actual rise and fall distributions shown in Figure [71 are suggestive of a broad, single-peaked distri- 
bution in RMF for the SDSS-II supernovae. The width of the RMF histogram is wide compared 
to the errors in the measured rise and fall times. This implies that the rise time is not strongly 
correlated with the fall. 




t, - tf (doys) t, - tf (doys) 

Fig. 6. — : Left: standard histogram of the rest-frame tr ~ tf distribution for 105 SDSS-II SNe. This figure is 
suggestive of only a single class of light curves, at least in terms of their tr — tf values. Right: the ideogram oitr~tf 
values, where each supernova is treated as a Gaussian and then these Gaussians are added together. This is also 
representative of a single distribution b\tr~tf, and agrees very well with the standard histogram. The ideogram is 
described in detail in Section 5.1. 

The ideogram in the right panel of Figure [6] is a visual device only and is used to display a 
histogram of the data weighted by errors m. — tf. It is not used to numerically evaluate any 
properties of the data. The distribution of a standard histogram, one where the data values are 
separated into rigid bins and then simply counted, can be very dependent upon the location and 
width of the bins. The ideogram, however, treats each data point as a Gaussian distributed about 
the mean value, so that the appearance of the distribution is only dependent upon the data values 
and their errors. For our ideogram, we weight the area that each RMF value contributes to the 
overall distribution by l/a. In other words, instead of normalizing the Gaussian distribution to 1, 
we normalize to l/a. For each supernova RMF value, this gives us the distribution: 

_{x-RMFif 



^Particle Data Group, [http://pdg.lbl.gov72007/reviews/| textrpp.pdf 
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Fig. 7. — : Left: the rise (dots) and fall times (crosses) vs. Ami5(_B) for the 105 SDSS-II supernovae. Fall times are 
directly related to Ami5(_B) and the crosses show the relation. These data are for B band only, and are shown to 
reinforce that averaging the stretches introduces no bias to the data. The typical error ellipse for the rise time or fall 
time is shown in the lower right. The ellipse results from the covariance of the rise and fall estimates, resulting from 
the error in measuring the time of maximum brightness. Right: here, the B- and 1^-band stretch estimates have been 
combined. Of special interest in this figure is the tendency for the slowest declining SNe to be among the fastest 
risers, however, there are very few slow-declining slow-rising SNe found in the sample. Also of note is a minimum 
rise time of about 13.5 days. 

N 

, where ai is the error in RMFi. The ideogram is then simply given by I{x) = Di{x). 

i=l 

The tr — tf values are not uniformly distributed with fall time (Ami5(i?) ). Figure [5] displays 
the RMF versus Ami5(i?) for the 105 supernovae. We see that there are few slowly declining light 
curves (Ami5(i3)< 1.0 mag) with large RMF values. In other words, many of the slowest declining 
SNe are among the fastest risers. This point is reiterated by the actual distribution of rise and fall 
shown in Figure [71 

For comparison with the SDSS-II data, we plot in Figure [8] the RMF trajectory of a 19.5-day 
rise-time fiducial curve with a single stretch applied to both rise and fall. This is currently the 
'standard model' of a SN la light curve and it does pass through much of the real data; however, 
the standard model cannot match the fast rising or very slow rising light curves in the SDSS-II 
sample. 
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Fig. 8. — : Plot of tr — tf against Ami5(_B), showing only the 105 SDSS-II supernovae. This figure displays the 
unexpected result that the slowest declining SNe tend to be the fastest risers. The dotted line shows tr — tf vs. 
Ami5(_B) relation if SNe la followed a single-stretch model with a rise time of 19.5 days and a fall of 15 days. The 
solid line shows the predicted t,. — tf for a series of Kasen light curve models that vary total radioactive nickel yield. 
The size of the points represents the error, so that larger points have smaller error. 



5.2. Average Rise Time and Decline-Rate-Corrected Curves 



Using our errors to weight tlie rise times of the sample, we obtain for the nearby set of SNe la 
an average rise time of 16.82 it 0.28 days (standard error) with a standard deviation of 1.77 days. 
For the SDSS-II SNe la, we obtain a weighted average rise time of 17.38 it 0.17 days (standard 
error) with a standard deviation of 1.8 days. These values are in general agreement, especially 
considering that the SDSS-II super novae and the l ocal sample have different selection biases. This 
result agrees very well with that of lStrovinkI (j2007l ). who found an average rise time of 17.44 it 0.39 
days in the nearby sample. 



In order to gain further insight into the SDSS-II light curves, it is beneficial to perform a 
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declin e-rate correction to the sample of light curves, similar to the procedure described by IStrovink 



(j2007l ). This correction is performed by dividing the entire time axis of each SN by the fall stretch 
found by the 2-stretch fitter. Based on the ideas of the single-stretch method, that every SN la 
light curve maintains the same relative shape between rise and fall, this approach should normalize 
the entire light curve to the shape of our template. 

We applied a decline-rate correction to the 105 SNe in our data set, by dividing the entire 
time axis by only, and found that the spread of the rise portion of the light curves increased. 
For the SDSS-II SNe la, the average decline-rate-corrected rise time is 17.32 it 0.12 with a standard 
deviation of 2.89 days. The increase in the spread of the rise times after applying the decline-rate 
correction is a departure from the expected result implied by a single-stretch method. 

This increase in the rise-time spread can be seen graphically in Figure [U which shows the result 
of applying the decline-rate correction to the 105 SNe la that pass our error cuts. This plot uses 
-B-band data only, along with B band only stretches, not the weighted average that was used in the 
RMF analysis. The SNe plotted in red have RMF < 2.0 days, while the ones plotted in blue have 
RMF > 2.0 days (the color selection for these supernovae is not meant to imply anything about 
the intrinsic colors). The plot on the left of Figure [9] contains the SNe with no stretch corrections. 
Even without stretch corrections, the data show the distinct light curve shape differences in SNe la 
that are also shown graphically in Figure El Without stretch correcting, we see clearly that most 
of the slowest declining supernovae are also the fastest rising supernovae; this effect becomes more 
apparent after decline-rate correction. The slow-declining, fast-rising light curves (in red) rise even 
faster, while the faster declining, slower rising light curves (in blue) rise more slowly. The dispersion 
in rise time increases from the uncorrected mean value, since the average rise times agree between 
the corrected and uncorrected distributions. This is more evidence that current SN la light curve 
fitting models are not able to acceptably characterize the full range of SN la light curve shapes. 
We have arbitrarily selected the dividing point at 2 days in tr — tf, which is the mean value, and 
we emphasize that the light curves show a continuum tr — tj values displaying minimal evidence 
for a two-group classification using this parameter. 



5.3. Comparing the 2-Stretch Fitter with Single Stretch 

In order to further examine the result that the 2-stretch fitter is a better representation of SN 
la light curve shapes, we selected a subset of 99 SDSS-II SNe la (from the 105 that passed our 
2-stretch only error cuts) that passed our error cuts using both the 2-stretch fitter and a single- 
stretch fitter, and calculated the value of the stretch-corrected light curve. In this analysis, we 
calculated over the entire region that was fit. For the 2-stretch fitter, this composite light curve 
has a value of 1140 with 951 dof. For the single-stretch fitter, this composite light curve has a 
value of 1515 with 1050 dof. Adding another parameter to a fitting procedure should be expected 
to reduce the value by approximately 1 (for each individual fit, so in this case 99) if that fitting 
parameter does not add any new useful information to the results. In our case, the 2-stretch fitter 
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Fig. 9. — : Left: the 105 SNe la used in this study, with no stretch corrections apphed and B-band data only. The red 
points represent SNe la with tr ~tf < 2.0 days, while the blue points represent SNe la with tr — tf > 2.0 days. Right: 
the same supernovae, with the entire time axis divided by the fall stretch as found by the 2 stretch fitter, showing 
the tendency for the slowest decliners to be the fastest risers. This analysis implies that the previously accepted 
correlation between rise and fall, e.g., a single-stretch model, is not supported by the data once the rise and fall times 
are disconnected in the fitting process. 



has a significant impact in reducing the value of the composite light curve, demonstrating that 
the 2-stretch fitter is a better model for representing SN la light curves than a single-stretch fitter. 



5.4. Fitting Simulated Light Curves 



As part of the SDSS-II supernova analysis, a sof tware package called SNANA was developed 
to analyze data and to create simulated light curves (jKessler et al.l l2009bl l . SNANA can create a 
set of supernova light curves with the same cadence, redshift distribution, and photometric signal- 
to-noise ratio as the SDSS-II survey, created from a variety of light curve models. In this section, 
we test the MLCS2k2 model and the single-stretch model against the real SDSS-II SN la data. 

We synthesized 2694 simulated light curves using a 17-day rise-time MLCS2K2 te mplate with a 



range o f A values and extinction that matches the analysis of SDSS-II SNe la done in lKessler et al 



(|2009al ) (see Section 6). These simulated curves have a range of z- values from 0.019 to 0.35 with 
an average of 0.19. We applied the same error cuts to these simulated curves as we applied to the 
SDSS-II data; the z-values after error cutting range from 0.019 to 0.24 with an average of 0.12. 
This simulated distribution after error cutting is shown in Figure [TOl It displays characteristics 
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that are present in the SDSS-II SNe as weh; it contains many SNe that are very slow decliners but 
are among the fastest risers. However, the simulation is lacking SNe that have relatively normal 
fall times yet are slow risers such as SN 1990N. 




Fig. 10. — : Left: distribution using SNANA simulated curves modeled with the MLCS2K2 17.0-day rise-time tem- 
plate. From an initial set of 2694 simulated light curves, 620 pass our error cuts, a very similar fraction as for the real 
data. Note that the MLCS2K2 model does not capture the true spread in rise times observed in the SDSS-II data. 
The rise time stays mostly constant, except in slowly declining SNe where it decreases with increasing fall time, an 
interesting correlation between rise and fall that is observed in the SDSS SNe as well. Right: simulated curves from 
the same input template but for a single-stretch model, and fit with the 2-stretch fitter. The rise time follows the fall 
time exactly for all values of decline rate, and there is no apparent trend for slow-declining light curves to have fast 
rise times in this single-stretch simulation. 



We also synthesized 2651 curves using a 17-day rise-time templa te simulated with a single- 
stretch model. The stretch values are based on the same analysis from lKessler et al.l (|2009al ). The 
z-values range from 0.019 to 0.35 with an average of 0.25. After error cutting, the z-values range 
from 0.019 to 0.34 with an average of 0.18. We fit these simulated curves with the 2-stretch fitter, 
but their underlying distribution is based on single stretch. The results can be seen in Figure [TOl As 
expected, this single-stretch simulation shows a strong correspondence between rise and fall that is 
not apparent in the SDSS-II data. Upon comparison with Figure [71 there is a noticeable difference 
in the spread of RMF values, with the SDSS-II data having a much greater range. Once again this 
implies a lack of correlation between rise and fall that is unachievable in this single-stretch model of 
light curve generation. However, the single-stretch simulation does produce some normal-declining 
events that have very slow rise times. 



Our simulated curves using an MLCS2k2 model and a single-stretch model imply that the 
range of SN la light curves is not fully matched with current light curve fitters. MLCS2K2 models 
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appear to be slightly better than single-stretch models as they capture some of the slow- declining 
fast-rising SNe la that are evident in the real data. Indeed we expect MLCS2K2 to be better 
simply because it is trained on real SNe. Also note that the MLCS2k2 synthesized curves maintain 
an approximately constant minimum for the rise time of fast-declining SNe at around 13-14 days, 
which appears in the SDSS-II SNe but could be a result of meager numbers of these SNe la. This 
could be an indication of a physical constraint on the physics of SN la explosions, as it seems that 
increasingly fast decliners do not correspond to increasingly fast risers. 



6. Discussion 



6.1. Comparison to Previous Rise-Time Studies 



We obtain an average rise time for the SDSS-II SNe la of 17.38 it .17 days (standard error 
of the mean); this is significantly different from the result obtained in iRiess et al.l (|l999l ). For 
the entire nearby sample, we obtain an average rise time of 16.82 it 0.28 days (standard error). 
The standard deviation in the rise time is 1.77 days. This is in general agreement with the value 
reported fo r the SDSS-II data. The reason that our results differ from the 19.5 days found by 
Riess et al.l (119991 ) is the difference in shape of the fiducial cur ves used in the two s tudies, and 
the application of the single-stretch fitter in previous methods. lAldering et al.l (|2000l ) prove that 
even slight variations in the declining portion of the template can change the measured rise time 
by 2 or more days using a single-stretch met hod. Figure [2] shows that the "Leibundgut" template 
(|Leibundgutlll989l ) used bv lRiess et al.l (jl999l ) differs greatly from the MLCS2K2 curves in the pre- 
maximum phase. The Leibundgut template is 0.66 mag fainter than maximum light at —10 days 
while the MLCS2k2 fiducial curve reaches 0.66 mag at —8.2 days. Such a broad template combined 
with a single-stretch fit tends to force the rise to be very slow and to push the e stimated time of 
maximum later than the true peak brightness (see Figure 7 in lRiess et al.l (|l999l )). The narrower 
MLCS2k2 fiducial combined with the 2-stretch fitter shows that the rise is faster than that found 
by several previous studies and that there is a wide range of observed rise times for a given fall 
time. 

We created extrapolations to the Leibundgut template in the exact same manner as for the 
MLCS2k2 template, ranging in rise time from 14 to 20 days in half-day intervals. Using the tem- 
plate determination method described in Section 4.3, we found the best extrapolated Leibundgut 
template had a rise time of 19.5 days (note that this utilizes the 2-stretch fitter). However, after 
our error cutting, we obtain an average rise time of 17.18 it 0.19 days using this template. Even 
though the extrapolated Leibundgut template has a longer rise time than the MLCS2k2 template, 
the stretches obtained are smaller so that the average is consistent. Indeed this should be expected 
with the 2 stretch fitter; the template with the best relationship between early rise and later rise is 
selected, but the overall rise time remains constant despite the difference in the fiducial relationship 
between rise and fall. 
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To fur ther demons t rate t his point, we also performed a rise-time extrapolation using the same 
method as iRiess et alj (|l999l ). as well as using the Leibundgut template. We fit all 391 SDSS-II 
SNe la with a single-stretch fit using only data from -10 days to 25 days, and we cut SNe with 
stretch errors greater than 2.0 days. After stretch correcting the full light curve, we used the data 
less than -10 days (which was not used in the fitting process) to fit a parabolic extrapolation to 
explosion. Using this method, we obtain an e xtrapolation of 19 .64 days for the SDSS-II SNe la. 
This agrees quite well with both the results of iRiess et al.l ()1999l ) and our own template selection 
using Leibundgut template extrapolations and the 2 stretch fitter. 

We have shown in this section that the difference between our rise time and that of previous 
studies revolves around the flexibility of the 2 stretch fitter regarding the template used. A single- 
stretch fitter that is used to determine the best extrapolation to explosion is much too reliant on 
the relationship between rise and fall in the template. The 2 stretch fitter is more flexible in the 
sense that the input template can have any relationship between rise and fall and the average rise 
time will be consistent with other templates. Performing a single-stretch fit and extrapolating to 
explosion gives the same result as our template determination method from Section 4.3. However, 
in order to find the average rise time, the addition of another parameter to independently estimate 
the rise stretch is required, because the rise and fall are not strictly correlated. 



6.2. Shape of the Early Rise 



As with previous studies, we have assumed that the optical flux rises as soon after explosion. 



ArnettI (119821 ) showed that adiabatic losses should nearly balance heating by radioactive decay and 
keep the effective temperature relatively constant. So the lui ninosity goes as the radius squared 
and therefore the time squared for a constant expansion rate. IConley et al.l (120061 ) directly tested 
this model by fitting light curves with the temporal power-law index as a free parameter. They 
found the temporal index n that best matched the low- and high-redshift supernovae was 1.8ib0.2, 
consistent with the Arnett calculation. 

To estimate the shape of the early rise, we fit the 105 SDSS-II light curves with good rise and 
fall data with the 2-stretch fitter using only epochs later than —10 days. We then corrected all 
the light curves to the fiducial curve using the rise and fall stretches which provides 103 i?-band 
observations between —20 and —10 days before maximum. We then calculate the parameter 
over this interval by comparing the data to the function / = A{t — to)"^, where the time of explosion 
(to) and power-law index (n) are allowed to float. The value of the parameter A is set by the 
condition that the function must meet the MLCS2k2 fiducial curve at —10 days. 



l.SOtoil- The index is 



We found the best-fit power-law index for the SDSS-II data was n 
strongly correlated with the time of explosion w hich we found to be -16.8+[j| days. The result 
is consistent with the Arnett prediction and the IConley et al.l (120061 ) estimate for the early light 
curve shape. 
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We can also test whether or not the early color of the SDSS-II supernovae follows Arnett's 
prediction of slow and modest temperature changes. For each SDSS-II observation where there is 
a good rest-frame B and V measurement, we construct a color based on the V — B normalized flux 
difference. This essentially compares the supernova color to the color at S-band maximum and 
avoids the problem of dealing with magnitudes at low flux levels. 

As a consistency check, for each supernova we take the flux color averaged between —15 and 
— 10 days to see how it varies with redshift. We find that the color is essentially constant as a 
function of redshift which gives us confidence that the SNANA /c-corrections at early times are 
self-consistent. We also see that the colors for z > 0.3 are very noisy and we restrict our color 
analysis to the 332 events at redshifts less than 0.3. 

Figure [TT] shows the flux color as a function of supernova phase. The colors have been binned by 
rest-frame day and the median calculated for each bin. As expected, the supernovae get significantly 
red after maximum as the T^-band light curve fades more slowly than the B. Near the time of 
explosion, the color is slightly red and shifts to the blue just before maximum. This observed color 
change is exactly what is expected given the 2-stretch fit of the individual B- and V-band light 
curves (and this is shown by the solid line in the figure) . Using the median i?-band fiux, we convert 
the fiux color into magnitudes and find the color index at —15 days is B — V = 0.5 which shifts 
to B — V = 0.0 at —9 days. The data earlier than —15 days is too noisy to attempt estimating a 
reliable magnitude colo r index. We convert the colo r index to an effective temperature using the 
empirical relation from Sekiguchi fc Fukugita (j200ol ) and find that T^ff rises linearly from 6000 K 



at —15 days to 9500 K at —9 days and then stays nearly constant until maximum. 



This linear rise in temperature over the first week after explosion is not expected from lArnett 



(j 19821 ). and is puzzling. If the color variation represents a temperature change and the optical 
bands are following the bolometric fiux, then we expect a temporal power-law index of at least 
n = 4 instead of the observed n w 2. In conclusion, we find that the power-law index for the SDSS 
SNe la agrees with the Arnett prediction; however, we find that the early light curves do not follow 
the prediction of slow and modest temperature changes. 



6.3. Source of Rise-Time Variation 



Our SDSS-II supernovae show a ran ge of rise time s for a fixed fall time, suggesti ng that the 
physics of the rise and fall epochs differ. lArnettI (j 19821 ) and iPinto fc EastmanI (j2000al ) show that 
the light curves of SNe la can be simply described by the deposition of energy from synthesized 
radioactive elements combined with the diffusion rate of energy out of the expanding nebula. Before 
maximum brightness the energy input rate from radioactivity exceeds the energy lost to luminosity. 
At maximum, the luminosity matches the instantaneous energy deposition rate and t he decline from 



maximum correlates with the total radioactive yield. Using more detailed models, IWoosley et al 



show that the shape of the light curve depends on more than the radioactive yield. For 
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Fig. 11. — : The V — B flux color relative to Umax as a function of SN phase. The data have been binned in 1-day 
(rest-frame) intervals and the median of each bin calculated. The diamonds represent the SDSS-Il SNe with redshift 
less than 0.3. The solid line represents the color of the template. The dotted line is the la spread of the data. 
Calculating a color index for times later than -15 days, we find that the effective temp erature i ncrea ses linearly from 
-15 to -9 days then remains nearly constant until max;imum. This is not expected from Arnett ( IQSd ). which predicts 
that the effective temperature should remain constant. This implies that the power-law index should be nearer to 
n « 4, rather than the n « 2 calculated from the observations. 



example, increasing the production of intermediate mass elements (e.g., silicon and calcium) at a 
fixed nickel mass (^^Ni) narrows the light curve while increasing the peak luminosity. This occurs 
because the total burned mass correlates with the KE of explosion, and the faster the supernova 
expands, the earlier the radioactive energy diffuses out. In contrast, a low KE means a slow rise 
and faint maximum. Note that the variation in KE trends in the opposite direction to the Phillip's 
relation. Observationally, the radioactive yield and KE influence the rise and fall times in different 
ways and variation in these two parameters could help explain the range of rise times seen in the 
SDSS-II supernovae. Other explosion parameters, suc h as non-radioactive iron yield or degree of 
mixing, could also impact the rise versus fall times (IWoosley et al.l l2009l ) . For illustration, we 
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consider only KE variation here. 



Kasen fc WooslevI (|2007l ) have calculated model light curves for SNe la with fixed KE and 
varying radioactive nickel masses to show that they can match the Phillip's relation fairly well. 
We can apply the 2-stretch method to these model curves almost as easily as to real data. From 
this we can see if the rise and fall times of the models can refine the physics of the observed 
light curve shapes. Unfortunately, the MLCS2k2 fiducial does not match the model light curves 
sufficiently well to estimate accurate rise and fall times. Instead, we chose the M^i = 0.49 Mq 
model curve as a fiducial and applied the 2-stretch fitter to the remaining models. The resulting 
rise times for fixed KE models with varying nickel yield are shown as a dashed line in Figure [12] and 
compared with observed low-redshift rise/fall times. The model rise time varies only modestly with 
Ami5(B) (slope of —2.6 days/Ami5(i?) ) compared to the fall time (slope of —9.2 days/Ami5(i?) 
at Ami5(i?) = l.l). This is in stark contrast to the single-stretch parameterization where the rise 
time parallels the fall time in this diagram. 

We compare the predicted effect on rise/fall times from varying ^^Ni yield with the SDSS-II 
data in Figure [HI The observed trend of tr — tf increasing with larger Ami5(i?) is well matched 
by the model, although there are still many supernovae in the sample that lie far from the model. 
Clearly, the single-stretch parameterization of light curve shape is not justified by the observed 
rise/fall-time variations or predictions of the ^^Ni model. The small scatter in the rise portion of 
the SN la light curves seen in the full SDSS-II sample (Figure [TJ should not be a s urprise if it is 
a dir ect consequence of ^^Ni yield being the dominant source of diversity in SNe la ( Arnett et al 
19851 ). This is a popular interpretation for variations in peak brightness in SNe la. 



Kasen &: WooslevI (|2007l ) models that vary KE at fixed nickel yield were also fit with the 2- 



stretch method. Arrows in Figure [T2] illustrate the amplitude and direction in the rise/fall time 
versus Ami5(i?) plane of a supernova that increases its nickel yield by 25% at fixed KE and then 
increases its KE by 25% at fixed radioactive yield. The change in KE produces a steeper variation 
in rise time than does nickel and in the opposite sense, i.e., increasing KE results in a faster light 
curve decay as well as a shorter rise time. Unlike nickel variations, changes in KE affect the rise 
time (slope of —9.7 days/Amis (i?)) as strongly as the fall time, meaning pure variations in KE 
result in almost no change in tr — tf (and would be well fit by a single-stretch parameterization). 

While this analysis focuses on the effects of KE and ^^Ni yield, there are other parameters 
that may be important to the light cur ve shape as well. In particular, mixing the ^^ Ni to a larger 
radius may result in fas ter rise times (iPinto &: Eastman! l2000bl : iHoflich et al.l 12010 ) . This effect 
plays a secondary role ( Hoflich et al. 201o|) in light curve shape, however, and is not included 
in this analysis. The iKasen &: WooslevI (|2007l ) models used in this analysis contain no effects 
from variations in ^^Ni distribution. We have only employed model light curves with fixed iron 
mass yields, which is th e parameter that controls the ^^Ni distribution in the one-dimensional 
Kasen fc WooslevI (|2007l ) models. 



In principle, varying both ^^Ni yield and KE will allow a supernova to reach any point on the 
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Fig. 12. — : Plot of rise time and fall time vs. Ami5(_B) for the 14 nearby SNe la analyzed in this study, with 
theoretical distributions derived from Kasen model light curves. The dotted line is a representation of the Kasen 
model rise times and its shallow slope implies that rise time is only mildly dependent on radioactive yield. The solid 
line is a quadratic fit to the fall times estimated from the models and are very consistent with the observed fall times. 
The arrows represent how a 25% change in KE or ^^Ni yield would change the position of a specific SN la from the 
Kasen models. 



tr — tf versus Ami5(i?) diagram (Figure E]). Changes in radioactive yield will move supernovae 
diagonally, while KE variations shift supernovae horizontally. Given these motions, it is difficult to 
explain the handful of events with very slow rise times and normal decline rates {t^ — tf > 6 and 
Ami5(i?)~ 1.0). SN 1990N may be a member of this group in the low-redshift set. The Kasen 
models also predict that peak luminosity should be rela ted to nickel yield ari d KE variations, with 
an increase in ^^Ni and KE resulting in brighter events. IWooslev et al.l ( 2003 ) found, however, that 
the effect of these parameters on light curve width works in the opposite sense. An increase in KE 
narrows the light curve while increasing peak luminosity and a larger nickel mass widens the light 
curve while making a brighter peak. Examining the absolute magnitude estimates of the SDSS-II 
events may help sort out the origin of the rise/fall variations. 
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6.4. Correlation with Host Galaxy Color 



Scannapieco fc BildstenI (j2005l ) and ISullivan et alj (j2006l ) have shown that rates of SNe la 
are strongly connected with host galaxy star formation rate. They model the rates as coming 
from two sources, one from a passive population of stars and the other a prompt population of 
supernovae correlated with high star formation rates. It is natural to speculate that the range 
of RMF corresponds to the sources of supernovae in the two stellar populations. To test this 
idea, we construct the g — i color of the host galaxies of our sample of 105 SDSS-II supernovae. 
The g — i color index correlates well with specific star formation rate measured from the SDSS 
spectroscopic sample. The host magnitudes are taken fro m the SDSS-DR5 cat alog and are color 
corrected to the rest frame by interpolating from tables in lFukugita et al.l (j 19961 ). The host galaxy 
color versus supernova tr — tf is shown in Figure [T3l for supernovae with Ami5(-B) < 1.5 mag. Very 



fast d e clining supernovae are known to be associated with early-type galaxies (e.g.. [Gallagher et al, 
(j2005l ): iHamuv et al.l (jl996al )) and might confuse the result. 



There is no significant correlation between host g — i color and the rise-time properties of the 
supernovae. Supernovae with RMF < 2.0 days tend to have slower decline rates than the other 
supernovae and can be expected to result from the prompt, high star-formation rate population. 
We do not observe this relation in our sample. 



6.5. Correlation with Hubble Residuals 

The Phillip's relation shows that the SN la light curve decline rate is related to the optical 
peak luminosity. The SDSS-II light curves suggest that there are a range of basic light curve shapes 
and this leads to several questions: 

• Which parameter is the best indicator of peak luminosity: fall time, rise time, total width, 
or rise minus fall time? 

• Is there a difference in the average peak luminosity between fast and slow risers for a fixed 
decline rate? 

• Do the color or dust properties vary with rise time? 

To investigate these questions, we estimated the apparent peak F-band magnitude, and peak 
color by fitting the 105 light curves with good rise/fall-time measurements with MLCS2k2. We 
relaxed all priors to the fits (fiatnegav) and assumed a fixed extinction law of Ry = 2.5. Despite 
the good rise-time information in the light curve data, we restricted the fit to include only those 
points later than —10 days before maximum; this was done because we found rather poor fits to the 
light curves when using the MLCS2k2 template extrapolations that put the explosion date at —20 
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Fig. 13. — : Host galaxy color vs. the difference in rise and fall times. Galaxies that hosted SNe la with t,- —tf < 2.0 
days are plotted as diamonds while hosts with high RMF SNe are circles. The size of the plotting symbol correlates 
with the total brightness of the galaxy. The RMF cut at 2 days was selected arbitrarily. The host color distributions 
do not appear different between these two groups. 



days before B maximum. For this analysis, we cut nine supernovae with Ami5(i?)> 1.5 because 
their colors are strongly dependent on decline rate (jGarnavich et al.ll2004l ). 



We also fit the light curves using the SALT-II software (|Guv et alj 120071 ) and compared the 
color parameter, c, and peak apparent magnitudes to those from MLCS2k2. The average color 
difference was {B — V) — c = —0.027 with an rms dispersion of 0.045 mag and the average apparent 
magnitude difference was my (mlcs)— my (salt) = 0.019 with an rms dispersion of 0.038 mag. We 
conclude that for these well-sampled light curves the differences in fit parameters between MLCS2k2 
and SALT-II are insignificant. 



Given the estimated apparent peak magnitude, redshift, {B — V) color, and our measured 
rise/fall stretches, we minimized the residuals on the Hubble diagram calculated from 
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Am = mv - 5log{dL) + 25 - a + f3{B - 1/) + 7(5/ - 1) + 6{sr - 1) + e(tr - */) 



where a, f3, 7, 6, and e are coefficients allowed to vary in the minimization. The luminosity 
distance {cIl in Mpc) is calculated from the known redshift and assuming the cosmological param- 
eters of Hq = 65 km s~^ Mpc~^, Q\ = 0.7, and = 0.3. Choosing a different set of cosmological 
parameters changes the relative absolute magnitudes by only a few percent at these redshifts. An 
'amoeba' algorithm was employed for the minimization of the parameter, 

i 

, where Sfii is the distance modulus uncertainty from the MLCS2k2 fit. We calculate for several 
combinations of the free parameters and the results of these fits are shown in Table [2j 

The scatter about the Hubble line for a one parameter fit (zero point) applied to the remaining 
96 events provides an rms scatter of 0.21 mag. This small scatter is not surprising given that 
the quality of the data is excellent and we have eliminated intrinsically red supernovae with the 
Ami5(5)> 1.5 cut. 

Adding a color term reduces the scatter to 0.18 mag. The size of the color term, /3 = —1.47, 
should correspond to Ry, the ratio between F-band extinction and reddening in B — V, if dust 
were the source of th e color variation. As noted for both local and high-redshift supernovae (e.g.. 



Kessler et al.l (|2009al )). the best-fit color term tends to be much smaller than the standard Milky 
Way extinction value of Ry = 3.1. This may indicate non-standard dust in other galaxies or that 
the supernova colors are not well understood. 

Minimizing the scatter by adding the fall stretch, sj, as a third parameter has some effect 
on the Hubble residuals, reducing the scatter to 0.16 mag. A fourth parameter that includes the 

Table 2. Hubble Residuals 



Fit 


96 SNe 


96 SNe 


96 SNe 


96 SNe 


96 SNe 


52 SNe (tr-t/ < 2.0) 


44 SNe (t,.-t/ > 2.0) 


Parameters 


1 


2 


3 


4 


3 


3 


3 


Mv (a) 


-19.087 


-19.093 


-19.067 


-19.066 


-19.061 


-19.054 


-19.070 


[B - V) 




-1.466 


-1.236 


-1.251 


-1.213 


-1.343 


-1.082 


Sf - 1 (7) 






0.614 


0.583 








Sr-1 (5) 








0.792 








Sr + Sf- 2.0 (e) 










0.691 


0.600 


0.765 


tj 


0.206 


0.175 


0.163 


0.145 


0.145 


0.143 


0.145 


xVdof 


214/95 


147/94 


131/93 


100/92 


101/93 


59/49 


42/41 



Note. — Residual = mv-51og(c(i)+25-a+,9(B-V)+7{sj-l)+5(sr-l)+e (tr+tj) 
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rise stretch, Sr, or the rise minus fall time, tr — tf, provides the smallest scatter of 0.145 mag. 
Surprisingly, the Hubble residuals are minimized when both the rise and fall stretches have large, 
positive coefficients. The Kasen models that vary both nickel yield and KE predict that increasing 
rise time (lowering KE) should reduce the supernova peak brightness, so the rise-time parameter 
would have the opposite sign from the fall-time parameter. Our hypothesis that KE has a strong 
influence on light curve shape and brightness does not match the observations well, but other 
variables should be explored. Varying the mass of stable iron elements appears to affect the rise 
time and brightness in the same directions as the KE variation so are also disfavored by this analysis. 

Applying a three-parameter minimization with the light curve characterized hy tr + tf results 
in the same scatter as a four-parameter fit with tr and as separate light curve shape indicators. 
This result suggests that it is the total width of the light curve that correlates with peak luminosity 
and not the rise or fall separately. Indeed, a three-parameter fit with tr — tj does no better in 
reducing the scatter than a two-parameter minimization with just zero point and color. 

When we divide the sample into a 'low RMF' group with tr — tf < 2.0 days and a 'high RMF' 
subset with tr — tj > 2.0 days there appears to be no substantial differences in the Hubble residuals. 
Applying a three-parameter minimization (zero point, color, and total stretch), we find the color 
term for the low and high RMF is only marginally different (see Figure [T4]) . We note, however, 
that the low RMF group has more blue supernovae than the high RMF group, with 9 of the 10 
bluest events having tr — tj < 2.0 days. The relation between luminosity and full width of the light 
curve also shows a little difference between the two rise-time divisions. 

The average absolute magnitudes of the two groups after color and light curve shape correction 
differ by less than 0.05 mag, implying that the presence of a range in rise times has very little direct 
impact on cosmological measurements. 



6.6. Estimating ^''Ni Yields 



Arnett's rule (jArnettI Il982l ) states that the bolometric luminosity of a SN la at maximum 
light is very close to the instantaneous energy being deposited by the synthesized radioactive e l- 
ements. This equality has been used by several studies ( e.g., IStritzinger LeibundgutI (j2005l ): 
Contardo. Leibundgut. Vaccal (|200d ): iHowell et al.l (|2009l )) to estimate the mass of ^^Ni created 
in the supernova. However, the time between explosion and maximum light is a key parameter in 
this calculation and it has often been assumed to be a constant 19.5 days. In this paper, we have 
found that the average rise time is shorter than that estimated in previous studies and it is not 
strictly linked to the decline speed. This will change the distribution of ^^Ni yields when compared 
wi th earlier assumptions about rise -time properties. We estimate ^^Ni yields using the prescription 
of Stritzinger Sz Leibundgut (j2005l ) but with a simplified technique for measuring the bolometric 
luminosity. We used the time from explosion to peak flux as our measured rise-time values. We 
have already calculated the extinction-corrected absolute l^-band magnitudes for the well-observed 
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Fig. 14. — : Absolute l/-band magnitude vs. peak color for fast-rising {tr — tf < 2.0 days; filled diamonds) and slow- 
rising {tr — tf > 2.0 days; open diamonds) SDSS-Il supernovae after correction for light curve shape. No difference 
in peak luminosity or color slope is apparent. Note, however, that the fast-rising events tend to be more blue at peak 
with 9 of the bluest 10 supernovae having tr — tf < 2.0. 



objects in the SDSS-II sample and we use these to scale an average type la spectrum which we in- 
tegrate over a wide wavelength range to get quasi-bolomet ric UVOIR l umino sity. For a supernova 
with absolute magnitude My =0.00 (Vega system) and the lHsiao et al.l (j2007l ) spectral template at 
S-band maximum, we find a bolometric flux of 2.04x10"^ erg cm~^ s~^ between 300 and 1000 nm. 
This approximation is best for normal type la and is not as accurate for sub-luminous events which 
tend to be more red than normal at maximum. However, IContardo. Leibundgut. fc Vaccal (120001 ) 
found that the V band is an excellent indicator of bolometric energy even at late phases when the 
spectrum has become red. Under these assumptions, the ^^Ni yield is simply 
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, where tr is the rise time in days and q is a constant describing the accuracy of Arnett's rule (we 
assume a = 1). 
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Fig. 15.—: Plot of the derived ^^Ni mass for the SDSS-II supernovae vs. the full width (rise plus fall times) of the 
B-band light curve. The ^^Ni mass is estimated from Arnett's rule so it is a function of both the peak bolometric 
luminosity and the rise time. The solid line is a least-squares fit to the points. Supernovae with slow rise times 
{tr > 18.7 days) are plotted as stars while those with slow declines {tf > 17.5 days) are shown as diamonds. Only 
one light curve had both a slow rise and decline and it is plotted as a triangle with a full width of 38 days. The 
supernovae with the largest ^^Ni yields (> 0.7Mq) have slowly rising light curves while slowly fading events produce 
more typical amounts ranging between 0.4 and 0.6 Mq of ^''Ni. Clearly, knowledge of the rise time is critical in 
understanding the physics of type la events. 



The estimat ed ^^Ni yields for the S DSS-II events range from 0.2 to 0.8 Mq, which are typical 
for this method ( Stritzinger et al.ll2006l ). Figure [15] shows the ^^Ni yield as a function of the total 
light curve width {tr + tf) and there is a strong correlation between these two quantities. Most 
of the supernovae producing the largest amount of ^^Ni are events with tr > 18.7 days (slowest 
rising 20%). In contrast, the slowest declining supernovae {tf > 17.5 days) produce average yields 
of ^^Ni. This demonstrates the importance of independently measuring the rise time of SN la to 



understand the physics of the thermonuclear explosion mechanisms. 
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7. Conclusion 

The SDSS-II supernova sample provides tight constraints on the rise time and shape of SN la 
light curves. The fc-corrected, median light curves of 391 events shows a i?-band rise time of 16.8 
days and a significantly smaller dispersion on the rise portion of the curve than on the fading side. 
This implies that th e rise time is less impacte d by variations in radioactive nickel yield than the fall, 



as predicted by the lKasen Sz Woosleyl (120071 ) model light curves. It is clear that the single-stretch 
parameter commonly used to characterize SN la light curves is not capable of capturing the full 
range of light curve shapes, especially when pre-maximum observations are available. 

We selected 105 SNe la from the larger sample that had sufficient photometric precision and 
cadence to yield rise- and fall-time errors of less than 2.0 days. Prom this set, we find the rise 
time in the B band is tr = 17.38 it 0.17 days (standard error of the mean ). This is signific antly 



different from the ri se time of 19.5 day s measured from local supernovae (|Riess et al.lll999l ) and 



high-redshift events (IConley et al.ll2006l ). We find the cause of this discrepancy to be the difference 
in the fiducial curves used in the analyses and our application of a 2-stretch fitting method that 
permits the rise and fall to be fit independently. We find rise times ranging from 13 days to 23 days 
for events with "normal" decline times of around 15 days, demonstrating that the rise and fall are 
not strictly correlated. 

We app l ied ou r 2-stretch light curve fitting method to data from SDSS-II to test the conjecture 



by lStrovinkI (|2007l ) of a bimodal distribution in SN la light curve shape. From our 105 high-quality 
light curves, we are not able to reproduce this bimodal distribution in the difference between rise 
and fall times. We do find a significant spread in RMF times for the sample, which implies again 
that the rise and fall are not strictly correlated. This large range in RMF is expected from explosion 
models that simply vary the radioactive nickel yield. 

Contrary to the premise of the single-stretch method, many of the slowest declining SNe 
are among the fastest risers. This effect is better modeled through MLCS2k2, as shown by our 
simulations, and this may be a result of MLCS2k2 being trained on a large number of real nearby 
SNe la. While MLCS2k2 is better than single stretch at representing the observed light curves, the 
simulations with MLCS2k2 do not capture the full spread in rise times evident for normal-declining 
SNe la. Specifically our simulations using MLCS2K2 do not contain normal-declining events that 
are slow rising. 

We investigated the relation between peak luminosity, light curve shape, and color by minimiz- 
ing the scatter on the Hubble diagram for the 96 SNe la with good rise/fall data and Ami5(i?)< 1.5 
. The Hubble residuals have the smallest scatter when using both rise and fall time as fit param- 
eters. The correlation of rise and fall time with luminosity have the same sign, meaning that the 
full width of the light curve is the best indicator of luminosity. Twelve models suggested that a 
second parameter, such as KE, competing with radioactive nickel yield, could be revealed in the 
rise and fall times and their correlation with luminosity. However, rise and fall times are found to 
act in concert in the SDSS-II SNe, implying that KE (or other physical parameters that broaden 



- 36 - 



the light curve and lower peak luminosity) is not a major contributor to light curve shape. 



Pignata et all (|2008l ) recently suggested that SNe la with high Si II velocity gradients (HVG) 
may have fast- rising light curves. We are analyzing the spectra of the SDSS sample to see if there 
is a rise-time, spectral velocity gradient correlation. It is interesting to note that SN 1990N, a 
low-redshift supernova with a very slow rise, has a low Si II velocity gradient, but a more complete 
sample needs to be constructed. 

The application of the 2-stretch fitting method we have developed provides a better fit to SN la 
light curves than simply stretching the entire time axis. By independently fitting the rise and fall 
portions of each light curve, we found that the rise and fall times are not strictly correlated. When 
correcting SN la luminosity for light curve shape, adding the rise-time parameter to the fall time 
reduces the scatter in the Hubble diagram when compared to using the fall time alone. However, 
we find that the full width of the light curve (rise plus fall time) minimizes the scatter about 
the Hubble line as well as using both rise and fall times as independent parameters. This does 
suggest some danger doing cosmology by combining a low-redshift supernova set dominated by 
poor pre-maximum sampling with a well-sampled high-redshift set. 

Our results show that measurement of rise time as a separate parameter characterizing SN la 
light curve shape may be a useful diagnostic in understanding the progenitors and explosion mecha- 



nisms of th ermonuclear events. Based on ^^N i yield calculations in lContardo. Leibundgut. k, Vacca 



(j200d ) and IStritzinger &: LeibundgutI (j2005l ). our average rise time would result in a 15%-20% re- 
duction in ^^Ni yield as compared to models that use a 19.5-day rise time. However, there are 
also some SNe la that have rise times longer than 19.5 days, resulting in a higher ^^Ni yield. It 
is clear that accurate nickel yield estimates require a good measurement of the rise time for indi- 
vidual events. Future supernova searches that provide well-sampled, high-quality data well before 
maximum light will reveal more about the nature of SN la. 
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